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Abstract 

Multiuser diversity (MUDiv) is one of the central concepts in multiuser (MU) systems. In particular, 
MUDiv allows for scheduling among users in order to eliminate the negative effects of unfavorable 
channel fading conditions of some users on the system performance. Scheduling, however, consumes 
energy (e.g., for making users' channel state information available to the scheduler). This extra usage 
of energy, which could potentially be used for data transmission, can be very wasteful, especially if 
the number of users is large. In this paper, we answer the question of how much MUDiv is required 
for energy limited MU systems. Focusing on uplink MU wireless systems, we develop MU scheduling 
algorithms which aim at maximizing the MUDiv gain. Toward this end, we introduce a new realistic 
energy model which accounts for scheduling energy and describes the distribution of the total energy 
between scheduling and data transmission stages. Using the fact that such energy distribution can be 
controlled by varying the number of active users, we optimize this number by either (i) minimizing the 
overall system bit error rate (BER) for a fixed total energy of all users in the system or (ii) minimizing 
the total energy of all users for fixed BER requirements. We find that for a fixed number of available 
users, the achievable MUDiv gain can be improved by activating only a subset of users. Using asymptotic 
analysis and numerical simulations, we show that our approach benefits from MUDiv gains higher than 
that achievable by generic greedy access algorithm, which is the optimal scheduling method for energy 
unlimited systems. 
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I. Introduction 

In wireless systems, unfavorable channel conditions remain the main hinderance to achieving 
desirable system throughput or bit error rate (BER). To overcome this problem in multiuser 
(MU) wireless systems, resource scheduling strategies, which use channel fading conditions 
as an opportunistic resource, have been proposed [QQ|, ED- Using the so-called opportunistic 
transmission, advanced scheduling strategies along with closed-loop designs 0, 01 have been 
developed in the literature (see J2|, 0, and references therein). The gain obtained by such 
opportunistic transmission methods is known as multiuser diversity (MUDiv) gain. 

For multi-point-to-point single-input single-output (SISO) wireless systems, the MUDiv gain 
was first studied in 0. The information theoretic results have shown that based on the optimal 
transmit power control, the overall system throughput can be maximized by allowing only the 
'best' user in a system to transmit at each time slot. For downlink MU systems, the MUDiv 
has been recognized as an effective method of improving the system performance measures 
such as spectral efficiency and quality of service over multipath fading channels [2]. As a result, 
MUDiv approaches have been adopted in commercial systems, e.g., systems based on orthogonal 
frequency division multiple access (OFDMA) 0. 

Toward improving the MUDiv gain, various system performance measures and their tradeoffs 
have been considered, as well as various algorithms have been developed 0- [fT3l . In |[8]|, the 
problem of multiuser downlink beamforming based on the MUDiv has been studied. In [9], the 
sum capacity caused by the MUDiv gain has been investigated with respect to two important MU 
system performance measures such as fairness and scheduling complexity. In [10], the delay- 
energy tradeoff in MUDiv systems has been analyzed. It has been shown that the energy required 
for guaranteeing an acceptable rate per user decreases at the cost of a longer delay. In [fTTI . |[T2~1 . 
low complexity scheduling strategies based on low rate channel feedback from users to the base 
station (BS) have been developed. In [fT3l . it has been argued that if the limited feedback is 
used, then the use of instantaneous channel norm feedback provides additional spatial channel 
information so that the MUDiv gain can be exploited efficiently in time, frequency, and space. 
However, in the existing literature, the MUDiv has been investigated for the case of fixed transmit 
resources, e.g., fixed transmit power and fixed number of active users. 

Although it has never been discussed before, it is important to note that the MUDiv gain 
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relies on the total energy available at the users and, therefore, depends on this energy, especially 
for the energy limited systems. Specifically, for scheduling purposes all users must share their 
own channel state information (CSI) with the BS per each data transmission. Then, a portion 
of the energy available at each user must be used primarily for scheduling, while only the 
remaining energy can be used for actual data transmission. Therefore, the important question is 
how to distribute the limited total energy available at the users between scheduling and actual 
transmission stages? If all users are active (available for scheduling) at all times, the waste of the 
energy used primarily for scheduling may be very significant. The latter will reduce the system 
performance. On the other hand, if only a small number of users is kept active per each data 
transmission, the corresponding MUDiv can be insufficient that also leads to system performance 
degradation. Therefore, the aforementioned question can be reformulated as the following signal 
processing question: how much MUDiv is required for MU systems? 

In this papeiQ, we develop methods which aim at maximizing the MUDiv gain in MU systems 
by exploiting a realistic energy model. Unlike existing schemes, we consider also the energy 
spent by users to make their CSIs available to the BS. By bringing this inherent energy usage 
into the picture, we find that it is better to choose (schedule for data transmission) the 'best' user 
from a subset of users (referred to as the set of active users) rather than among the entire set of 
users. The intuition is that, if a small subset of users is required to send their CSIs to the BS, 
more energy can be saved for actual data transmission and better overall system performance can 
be achieved. This is especially true for the energy limited systems. Thus, there is an inevitable 
tradeoff between the MUDiv and the energy saved for actual data transmission. Using this 
tradeoff, we aim at finding the optimal size of the set of active users so that either the total 
system bit error rate (BER) or the total energy of the users is minimized under practical system 
constraints. Using asymptotic analysis, we also study how much MUDiv can be achievable in 
various special cases of interest. 

The paper is organized as follows. The system model is introduced and the problem is described 
formally in Section EL In Section III, the systems performance measures such as BER, upper 
bound on BER, and approximate BER are derived. Section IV contains the answer to the main 
question of the paper, that is, how much MUDiv is required for MU systems, while Section V 

'Some preliminary results of this work have been published in 1151 . 
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provides some further analytical analysis. Extension to the case of multiple antenna MU systems 
is given in Section VI. Section VII presents numerical results and is followed by conclusions in 
Section VIII. 

II. System model and problem description 

A. System model 

Let K mobile users communicate with the BS. It is assumed for simplicity that each user 
as well as the BS is equipped with a single antenna. Thus, we consider an MU SISO system. 
This assumption, however, will be generalized to the case of multiple antenna MU systems in 
Section |VI] and it will be shown that such generalization is straightforward. 

Suppose that the wireless channel between user k and the BS is flat fading. The received 
signal at the BS from user k can be then represented as 

Xk = h k s k + v k , k = l,---,K (1) 

where the information-bearing symbol s k is a Gray-coded quadrature amplitude modulated 
(QAM) symbol] from a fixed constellation of size M, v k is the complex-valued zero-mean 
additive white Gaussian noise (AWGN) with unit variance, i.e., v k ~ CjV(0, 1), and h k ~ 
CJ\f(0, a k ) is the channel gain between user k and the BS. We assume that h k , Vfc are independent 
and known perfectly at the BS. 

One of the main concerns for scheduling in the heterogeneous MU environments is fairness 
among users. Among various fairness notions such as, for example, average throughput per 
user 0, variance of short-term throughput per user O, users' channel accessing period 031, 
our concern, in this paper, is fairness in terms of the equal user's probability of accessing the 
channel. According to this fairness notion, the scheduling is called fair if the channel accessing 
probabilities are equal for all users in the MU system. To satisfy such fairness conditions, we use 
an opportunistic scheduling (OS) scheme proposed in JTJ. This scheme incorporates an average 
power control which is instrumental for our further considerations of the energy distribution 
between scheduling and transmission stages. According to this scheme, a ratio of the actual 
signal-to-noise ratio (SNR) to its own average is used for both scheduling and data transmission. 

2 Note that the approach can be easily extended to other modulations. 
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The aforementioned scheduling scheme (hereafter referred to as greedy access (GA) scheme) 
gives equal chance to all users for accessing the channel. Thus, we employ it in this work. 

Let user k employ the average power control of [1J assuming that the variance of h k , i.e., o\, 
is known to him. Thus, the power is allocated to symbol s k in (QQ) so to obtain a desired average 
receive power at the receiver which must be the same for all users. Then, denoting the desired 
average receive power for unit transmit power by to, the corresponding transmit power at user k 
can be written as 

A fc = 4 A (2) 

where A is the transmit power before employing the average power control and tu\ is the 
desired average receive power that is equalized for all users via the average power control 
ojja\. Therefore, using © and instantaneous channel gain |/ifc| 2 , the instantaneous receive SNR 
at the BS from user k can be written as 

Pk = \h k \ 2 \k- (3) 
Since the variance of the AWGN in (OQ) is unit, © can be equivalently written as 

p k = \h k \ 2 uj\ (4) 

where h k ~ CJ\f(0, 1). Therefore, VA;, the distribution of p k is the same. 

Using © as a scheduling metric, we consider the GA scheme, where at a given time slot, 
the BS chooses only one out of multiple users for transmission. The user selection criterion is 
based on finding the user with the most favorable channel gain versus its own average. That is, 
user k* is scheduled for data transmission if 

k* = argmaxpfc- (5) 

k 

In practical MU environments, the system resources such as the number of users in the system 
and the total energy available at each user are usually limited. Under such system limitations, 
one interesting question is how the existing limits on the total available energy of all users should 
change the requirements on the MUDiv of the system. Indeed, one of the well known access 
schemes, i.e., the random access (RA) scheme (see lfT9lO . suggests to select users for transmission 
randomly one at a time. This scheme provides no MUDiv, i.e., K = 1, and, therefore, requires no 
extra energy spending for extra communications between the users and the BS at the scheduling 
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stage. On the other hand, the GA scheme improves the performance of MU systems due to its 
ability to select the 'best' user for transmission from the entire set of available users of size K 
HI, 0. The MUDiv of the GA scheme is then K = K. Unfortunately, in this case, the BS has to 
know the CSIs of all users in the system in the scheduling stage which requires additional energy 
spending. Therefore, if the total energy of all users in the system is limited, the use of the GA 
scheme may be very wasteful in terms of the energy spent at the scheduling stage. It reduces the 
energy available for actual data transmission that can lead to the system performance degradation. 
Therefore, the main query of this work is how much MUDiv K is required to improve the MU 
system performance? In other words, how many users should transmit their pilot symbols so to 
make their CSIs available at the BS. Based on these CSIs at the BS, one of the users is selected 
to access the channel. 

The aforementioned query can be solved by finding an optimal energy distribution between 
scheduling and data transmission, i.e., by selecting the cardinality K = \A\ G [1, K\ of a subset 
of active users A which participate in the scheduling. Here | ■ | denotes the cardinality of a 
set and the elements (users) of A are selected randomly in the beginning of every time slot 
according to a uniform distribution. Such random selection in each time slot is considered in 
order to achieve fairness among users in terms of equal channel accessing probability. 

Toward this end, let us first write the energies used for scheduling and data transmission as 
functions of K. Taking into account the scheduling stage, the energy consumed by user k at 
each time slot for both scheduling and data transmission can be defined as 

E T ,k = E 8 , k + l(k = k m )E d , ki Vfc (6) 

where E Sik denotes the energy spent for scheduling, E djk = T s X k is the energy spent for data 
transmission, T s stands for the symbol duration, and l(k = k*) is the indicator function which 
is equal to 1 if k — k* and otherwise^ Then, the total energy of all users over the time 
interval during which the average channel gains remain constant can be found as the sum of 
Ex,k, Vfc over many time slots N covering the whole interval. Since all users have equal chance 
of accessing the channel, at a given time slot, any user has access to the channel with probability 
1/ K. Then, it can be found that during N time slots, the energy used by each user for individual 

3 Note that without loss of generality, E s> k = E s j, Vfc, j is assumed. It corresponds, for example, to the practical situation 
when the codeword length of the transmitted symbol is long, while the number of pilot bits is relatively small. 
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data transmission is AfcT s • N/K where is the transmit power at user k which is equal to 
the kth user SNR under the assumption of the unit variance of the AWGN in (Q]). Therefore, 
asymptotically for large N, we can write that for K < K, the total energy is 

E K-GA A N (AiTs + . . . + XrTs) /R + KEf = Ed + REf (?) 

where the superscript ( y -) K ~ GA stands for the GA among K active users, Ef = NE Sjk denotes 
the energy consumed by each user for scheduling over N time slots, KEf denotes the total 
energy consumed by all users for scheduling, and Ed stands for the energy used by all users for 
actual data transmission. 

Although © is an asymptotic result, it is applicable to practical setups. Consider the random 
variable X corresponding to the actual number of time slots that a user is accessing the channel 
over iV time slots. Then X has a binomial distribution with average nix = N/K and standard 
deviation <r x = y/N/K(l — T/K) « y/N/K for large K. Therefore, for a x < mx/10, we 
need, N/K > 100 which is the realistic case in practical setups. 

B. Problem description 

Two different objectives can be considered for selecting K: (i) minimization of the system 
BER and (ii) minimization of the total energy consumed by all users in the system. Although the 
users are not connected to the same energy source, given the finite energies at individual users, 
the sum of individual user energies also determines the total energy consumed by all users. It 
is worth stressing here that for system performance analysis in MU systems, the total energy 
consumed by all users is more important than individual user energies because the MUDiv gain 
depends on the number of users participating in scheduling, and the energy which determines the 
MUDiv gain is the total energy consumed by all users, rather than the individual user energies. 
In addition, assume that for given channel statistics of all users, the energy consumption by 
each user over given time slot(s) is fixed on average. Then, the individual user energies are 
also fixed fractions of the total energy of all users on average (see JT), ED, [0, and references 
therein for similar observations for power or data rate). Since the individual user energies are 
fixed fractions of the total energy, by minimizing the total energy, the individual user energies 
are also minimized. 
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1) System BER minimization: In this case, we aim at minimizing the system BER under a 
constraint on E% ~ GA . Therefore, E% ~ GA is a constant independent of K, and it is straightforward 
to see that E d in © as well as KEf depend on K since the energy distribution between E d 
and KEf must be optimized by selecting such K that minimizes the system BER. Thus, for 
a given total energy consumed by all users, we first express the tradeoff between E d (K) and 
KEf as a function K. Let us define the ratio a = E GA /Ef where E GA denotes the energy for 
data transmission consumed in the generic GA scheme that holds K = K during all time slots. 
Here, the superscript (-) GA stands for the GA scheme. Then, representing Ef in terms of E GA , 
E^~ GA can be expressed as 

E«~ GA = Ka~ l E GA + E d . (8) 

Due to the fact that in the generic GA scheme K = K in all time slots, the total energy consumed 
by all users is a constant (denoted by E GA ). Constraining ([8]) to be equal to E GA , we obtain 
under such energy constraint 

KoT l E GA +E d = Ka~ l E GA +Ef A (9) 

--KE f =KE f 



pK — GA 



--E% A 



where the two terms on the right hand side represent the total energy E GA consumed in the 
generic GA scheme. It can be seen from © that if K is selected such that K < K, then more 
energy remains after scheduling, i.e., E d (K) = E GA — KEf. This extra energy can be assigned 
for actual data transmission, and E d (K) can be expressed in terms of K as 

E d (K) = ((K - K)a- 1 + 1) E GA . (10) 

Therefore, E d (K) benefits from the energy gain of (K — K)oT x + 1 if K < K. On the other 
hand, assigning more energy for scheduling KEf increases the MUDiv gain. Therefore, there 
exists a tradeoff between E d (K) and KEf, and the question now is where to spend the available 
energy in order to minimize the total system BER. One of the possibilities is to find the optimal 
value of K < K which minimizes the total system BER, while satisfying the constraint on the 
limited total energy of all users. 

2 ) Minimization of the total energy consumed by all users in the system: In this case, we 
aim at minimizing the total energy consumed by all users under the constraint that the system 



April 16, 2010 



DRAFT 



9 



BER remains below a pre-determined threshold. In order to satisfy the BER requirement, E d in 
CO) must remain constant (i.e., E d = E GA ) for any number of active users K. Therefore, the 
total energy E^~ GA can be minimized by selecting such K which minimizes KEf. Therefore, 
in this case, the total energy E^~ GA is also a function of K. More precisely, since E d = E GA 
for all K, then, by selecting K, E^~ GA can benefit from saving the energy at the scheduling 
stage. Therefore, E^~ GA can be expressed versus E GA as 

E*- GA = KE f + Ef A = (Ka~ l + 1)E$ A . (11) 

Since E GA = (Ka' 1 + 1)E GA and E GA is a constant in the considered energy minimization- 
based problem, E GA in CCD) can be expressed via E GA as E GA = (Ka' 1 + 1)~ X E GA . Using 
this relationship and (ITT]) . E^~ G can be further expressed versus E GA as 

E K-GA = Ka ~ 1 + 1 E GA_ (12) 

It can be seen from (fT2l) that the energy saving gain is (Ka^ 1 + l)/(Ka~ l + 1) if K < K. 
Therefore, the smallest possible subset of available users which satisfies the system target BER 
requirements is optimal in terms of providing the minimum E^ G . 

In order to express the aforementioned problems of selecting optimum number of active users 
formally, we first need to find an expression for the system BER as a function of K. 

III. System performance measures 

In this section, we derive expressions for the exact, upper bound (UB), and approximate 
BERs. The exact BER expression provides the highest accuracy for choosing K. However, it 
may require intense computations, which may not be practical in real-time. Therefore, a simple 
UB expression for BER is derived. The use of the UB BER expression instead of the exact 
BER in our problem will guarantee that the system BER requirements will be satisfied, but the 
resulting K may be sub-optimal. Therefore, approximate BER expressions, which require the 
minimum computations, are also derived. 

A. Exact BER expression 

The exact BER of the M-ary modulation over the AWGN channel can be written as [fT6l 

Pv b (M, p) = J2 CmMV^p) ( 13 ) 

1=1 
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where Q(x) = n 1 J^ 2 e li / 2sm2fl d^ is the error function. For a Gray-coded square M-ary 
quadrature amplitude modulation (QAM), the constants Q M , C M ,i, and c M ,i can be found in 

For a given to, the average BER is given by 

poo 

Pr 6 = / Pr b (M,y)p p (y)dy (14) 
./o 

where p p (-) is the probability density function (pdf) of p. 

Let p K ~ GA denotes p for K-GA scheme, i.e., p K ~ GA = p fe *. Considering the average power 
control, pk, Vk are independent and identically distributed (i.i.d.) random variables. Using this 
fact and applying higher order statistics, the pdf p p K-GA(y) can be found, for a given K, as 

e ~y/ n / / y\ \ K-l 



P^W-K—l^l,^)) (15) 
where = tu\, and 7(1, x) = (1 — e~ x ). Using (fl"5l) . the average BER can be written as ifTTl 

©A/ /•§ /"OO 

Pr«- GA (K) = J2C M ^- 1 K a e"* ( 7 (1, at))*" 1 dtd0 (16) 
i=1 -/o ./o 



where a = (Qgg + 1) and gg = CM,i/2 sin 9. Moreover, using the expression (3.312.1) in 11181 
p. 305], after some algebraic manipulations, we obtain the following closed form expression for 
G6J: 

e M 

Prf ~ GA (K) = V C M ^- l K / B(K, 1 + fl fi)d0 (17) 
i=i Jo 

where B(x,y) = Lt x ~ x {l — i) y_1 di denotes the beta function. For a given oj, it is clear from 
(fT7T) that Pr^ e ~ Gj4 (i ; r) depends on both .K" and A. Note that B(x,y) decreases exponentially with 
x at a given y. Therefore, for a given A, the system BER in (flTT) decreases exponentially with 
respect to K due to improvements in the MUDiv at the cost of increased E^~ GA in ©. 

B. Upper bound expression on BER 

The finite range of the integral in (fl6l) can be eliminated by considering the minimum value 
of gg. Thus, substituting 9 = n/2 in (fl6l) . we find the UB expression for (fT7T) with g u = Cjvf^/2 

4 Note that the BER of a Gray-coded coherent M-ary phase-shift keying (PSK) modulation in AWGN channel can also be 
expressed using l !13t . However, for brevity, only M-ary QAM modulation is considered here. 
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(Chernoff bound) as 

e M 

Pr«- GA (K) < Cm^K B(K, 1 + g u Q) = Pr«~ GA (K). (18) 
i=i 

This UB BER is clearly less complex than (fTTT) since it does not contain integration. 

C. Approximate BER expression 

Inserting (fl"5l) and the following approximation of (fl"3l) (6]|: Pr&(M, p) m 0.2e~ 9ap , where 
g a = 1.5/(M — 1), into (PT41) . the approximate BER expression can be written as 

Pr«- GA (K) « 0.2 X B(tf, 1 + <7 0) = Pr£- GA (if). (19) 

Note that in comparison to the exact and UB BER expressions, which have multiple summation 
terms of beta functions, the expression (IT9T ) requires minimum computations with a single beta 
function. 

Fig. CD illustrates the exact, UB, and approximate BER's, i.e., (fTTT) . (fl8l) . and (fl9"l) , of the i^-GA 
scheme for K = K = 1, 10, 50. It can be seen from this figure that the UB and approximate 
expressions produce the BER curves which lay within 0.5 dB of the exact BER. 

IV. Optimal selection of the number of active users 

Two scenarios are considered in this section for selecting the number of active users K 
optimally: (i) minimizing Prjf~ GA (-) while E^~ GA remains constant and (ii) minimizing E^~ GA 
while Pr^~ GA (-) is constrained to be acceptably smallO Optimization problems for each scenario 
are provided. 

The set of candidate values of K is the set of all positive integers smaller than or equal to 
K, i.e, /C = {1, • • • , K}. Note that due to hardware design limitations /C can be just a set of 
some integers smaller than or equal to K. The latter case can be easily adopted in the methods 
developed further. 

5 The exact, UB or approximate BER's can be considered. We use the notation ¥r^~ GA (■) to refer to any of these three BER 
expressions, i.e., Prf" G - 4 (.) G {Pr«-° A (K), Pr^- aA (K), V^~ GA {K)}. 
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A. Optimal selection of K based on system BER minimization 

Using © and (flOl) . under the finite energy constraint, we can find that the achievable energy 
gain for E d (K) determines X(K) in Prf ~ GA (-) as follow^ 

X(K) = E d (K) c/(NT s ) = ((K - K)a~ l + l) E G d A c/(NT s ) (20) 

where c = K(%2k=i a k 2 )~7 W - 11 is worth mentioning that E^ A /(NT S ) in d20]) stands for the 
average transmit power over N slots when K = K. Thus, it can be denoted as \ GA . Using this 
notation, (|20l) can be represented as 

\{K) = {(K-K)a- l + l) c X GA . (21) 

It can be seen from (|2TI) that a power gain of [(K — K)a^ + l) c in Prf" GA (-) is achieved 
if K < K. Denoting Pt = E/(NT S ) as the total average power consumed by all users during 
one slot, \(K) can be also expressed in terms of Pt as X(K) = c(Pt — -ft'a~ 1 A Gj4 ). 

Using \ (K) in (|2T|) . we optimize K for a given tu to minimize the system BER while satisfying 
the finite energy constraint. The corresponding optimization problem can be mathematically 
formulated as 

K* b = arg mm Prf ~ GA (iT) subject to E^~ GA = E. (22) 



One way to solve (|22|) is to employ a binary search over K e JC through direct computation 
of Prf ~ GA (-). However, direct evaluation of Prf ~ GA {-) is computationally complex, and such 
an approach can be inaccessible for applications sensitive to high computational complexity. 
Therefore, an approach, which avoids direct computation of Prf ~ GA {-) for all K, is proposed. 

To this end, let us relax K to be a real numbeiQ such that K G Let us also define 

r]{-) = T^Prf ~ GA (-). It can be observed that Prf ~ GA (-) is convex with respect to K due to the 
fact that J^Prf- GA (-) > 0. Thus, the minimum of Prf~ GA (-) over K 6 [1,K] can be found 
by minimizing f° r a given normalized fi, i.e., £In = wA GA . Note that this minimum is 

unique. Therefore, denoting K* as a real-valued solution, the corresponding optimal solution 
can be given by 



K* = arg min 

Ke[K n ,K!} 



V (K) 



(23) 



6 The argument K is added here to emphasize that A is a function of K in the BER minimization-based problem. 

7 While relaxing K to be a real number, we also assume that Pr^ _GA (-) is continuous on K. and differentiable at all points 



on /C. 
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Recall that /C is a finite set of integers, and the optimal Kl G /C may not be equivalent to K* 
in (|23T) . Therefore, (T23l) should be reformulated as 



K* b 



K, ifsign(7?(7r = K)) = -1 

1, else if sign(r/(iT = 1)) = 1 (24) 

^argmin^ e{A - ijXi+l} Prf ~ GA (K), otherwise 

where sign(a) = \a\/a for a G R with sign(O) = 1, and Ki G /C is the largest integer smaller 
than if* that satisfies the equality sign(?7(_K"j))sign(?7(i^ i+ i)) = —1. 

In order to find Kl G /C, we first compute sign^i^)) at K = K (and/or K = 1). If the 
resulting sign(-) is —1 (or 1), then we select Kl = K (or 1). Otherwise, Ki G /C can be found 
by binary search algorithm followed by selecting Kl at whichever of Ki or K i+ i that has a 



smaller Prf" GA (-). 



Considering, for example, the case when Pi b (■) = Pr 6a (K), it is shown in Appendix 
that 

1 f{K) 1 1 

^^-EzTT-ZWEiTT^yTT (25) 

where /(X) = g a uj\. 

Given VL N and a finite set JC — {1, 2, • • • , r/(i^) is illustrated in Fig. [21 It is shown that 
neither K — > 1 nor K K may minimize Prf : ~ Gj4 (-). Therefore, for a given f2jv, there exists 
an optimal 1 < Kl < K minimizing Pijf~ GA (-). For example, it can be seen from the figure 
that when M = 4 and K = 100, Kl rj 67 minimizes Prf ~ GA (-) at Vl N = A dB. 

5. Optimal selection of K based on E^~ GA minimization 

If an MU system is capable of recovering properly the transmitted information as long as the 
system BER is less than or equal to a predefined desired level, the optimal K can be found 
via minimization of E^~ GA in © subject to the constraint Prf ~ GA (-) < BER t where BER t 
is the required target BER. In this problem, different from the previous problem, E^~ GA is the 
optimization variable, while A (equivalently Vt^ = fi) is fixed. 

Two cases of delay tolerant and delay sensitive systems are of interest. 
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1) Delay tolerant (DT) systems: The constrained optimization problem for finding optimal K 
can be written in this case as 

K* t = arg min E^~ GA subject to Pv^~ GA (K) < BER t . (26) 

KaK. 

Note that for a given Vt, Pr^~ GA (-) and E^~ GA are, respectively, monotonically decreasing and 
monotonically increasing functions of K. Thus, among all values of K satisfying the constraint 
Prf ~ GA (-) < BER t , the smallest K e JC which minimizes E%~ GA is the solution of (|26l 

In order to find K dt G /C, we first need to find Pijf~ GA (-) when K = K. If K does not satisfy 
the system BER requirements, then K dt = 0. When K dt = 0, the system may allow delays to 
prevent the waste of the total energy of all users. Otherwise, the smallest K < K, which satisfies 
the system BER requirements, can be searched efficiently using, for example, a binary search 
algorithm. 

2) Delay sensitive (DS) systems: DS systems allow to transmit data even if Pr^~ GA (-) > BER t 
when K = K. Then, the corresponding constrained optimization problem can be written as 

argmin^ eK ; E^~ GA subject to Pr£~ GA (K) < BER t , if Prf~ GA (K) <BER t 
K* ds = { (27) 

K\, otherwise. 



The problem (1271) can be solved similar to the previous one. The only difference is that K ds = K 
even if K = K is not sufficient to satisfy the constraint Pr^~ GA (-) < BER t . 

V. Asymptotic analysis 
A. Asymptotic analysis of optimal K based on E^~ GA minimization 

In general, the optimum K based on E^~ GA minimization under fixed system BER cannot be 
found in closed form. However, its asymptotic behavior can be studied analytically. Recall that 
the derived system BER expressions depend on the beta function B(-, •). Thus, we first study 
the asymptotic behavior of B(-, ■) with respect to K. 

The following theorem summarizes the asymptotic behavior of the beta function. 

Theorem 1: Let x and y be two positive integers. When x — > oo, we have 

lim x y B(x,y) = T(y) (28) 

X'— >oo 

where T(y) = / °° t y-1 e~*di denotes the complete gamma function. 
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Proof: The beta function can be alternatively represented in terms of the following ratio of 
complete gamma functions 

B(«,») = 3?2M. (29) 
T(x + y) 

Using the expression T(x) — (x — 1)! for the gamma function, we can find the following ratio 

r(s) = (g-i)i = i (30) 

F(x + y) (x + y — 1)! (x + y — l)(x + y — 2) • • • (x + l)x 
Since (x + y — l)(x + y — 2) ■ • • (x + l)x in (l30l) is dominated by the first power term x y , the 
ratio in (|30l ), for x — )• oo, becomes 

lim r(x)/r(x + y) = x" 2 '. (31) 

x— >oo 

Thus, when x — >• oo, inserting (|3TT) into (|29l ) reveals the asymptotic behavior of (|29l ) as 

lim B(x, y) = lim T(y) = x~v F{y). (32) 

Since for a given y, T(y) in (|32l) is fixed, (|28T) is obtained when x — )■ oo. ■ 
Theorem Q] enables us to evaluate the system BER for two asymptotic cases of (i) large K 
and (ii) high SNR. We also aim at investigating how the optimal K scales asymptotically. For 
simplicity, only the approximate BER expression (fl9l) is considered in the further analysis. 

For the case of large K, we aim at analyzing Prjf~ GA (K) versus K and SNR. As per Theorem 
[H for large values of x, the following approximation holds true B(x,y) ~ x~ y T(y). Then, when 
interpreting x and y in (1281) as K and 1 + g a Vt, respectively, the system BER can be expressed 
for large values of K as 

Pr* - GA (K) = K-^ Q 0.2 T(l + g a Q) 

(33) 

= (K~ SNR ) 

where we use the alternative notation SNR= in the last expression. Therefore, (1331) shows 
how the system BER scales with respect to the MUDiv gain if K is large. 

For another asymptotic case of large SNR, Pt^ ga (K) can also be expressed in terms of 
SNR and K. Specifically, using the fact that B(x, y) = B(y, x), it follows straightforwardly from 
Theorem Q] that for large y, B(x,y) ~ y~ x T(x). Therefore, the system BER Prf ~ GA (K) can be 
expressed for large SNR as 

P^- GA (K) = n- K 0.2g- K T(K+l) 

(34) 

= 6 (SNR"*) . 
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It follows from (|34l ) that the total system BER scales inversely with the order of the SNR, i.e, 
the MUDiv is equal to K. 

Using (l33l and (|34l) . we can find the optimal K* E /C, i.e., either K dt for the delay tolerant 
or K ds for the delay sensitive systems. 

In the case when K is large, it can be found using (|33l) that for given to, A, and BER t , the 
optimal K* e /C, i.e., K dt for the delay tolerant or K ds for the delay sensitive systems, must 
satisfy the following inequality 

K* > (0.2r(l + g a ty)^ BER^ 



6(BER- 1/SNR ^ 



(35) 

It follows from (1351) that for large K and a given SNR, K* is an exponentially decreasing 
function of BER t . 

In the case of high SNR, it can be found from (|34l) that the optimal K* , which guarantees 
that the target BER is archived, i.e., the constraint Fxf~ (K) < BER t is satisfied, must obey 
the following inequality 

K , log BER,- = /logBER^X 
- logg„fi V log SNR I 



For a given BER t , it follows from (1361) that the corresponding optimal K* , i.e., K dt for the 
delay tolerant or K ds for the delay sensitive systems, is proportional to the inverse of log SNR. 
Moreover, unlike the case of large K, in the case of high SNR, K* decreases in a log-scale with 
BER t . 

B. Asymptotic analysis of optimal K based on the system BER minimization 

We again consider two cases of (i) large K and (ii) high SNR and study the asymptotic behavior 
of optimal K%, i.e., we study the asymptotic behavior of the solution of the optimization problem 
(|34|) . For simplicity, but without any loss of generality, we assume that c = 1. 

In the case when K is large, we first determine how the system BER scales with K while 
satisfying the finite energy constraint. The corresponding power gain given by (|2TT) is G p = 
((K — K)^ 1 + l) . Using this notation and (l33l) . the system BER can be asymptotically ex- 
pressed as 

Pr* - GA (K) = 0.2r(l + g a G p n N )K-*> G * n » 

(37) 

= e (k-^ snr ) 
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where we use the alternative notation SNR= Cl N . Therefore, if K ^> K — K, the achievable 
MUDiv gain is determined by G P (K) ■ SNR = ((K - K)a' x + l) ■ SNR instead of SNR. It is 
also worth mentioning that, for a given SNR, the asymptotic system BER scales exponentially 
with K^ R - K ^ a . The latter means, in particular, that the achievable system BER is lower in the 
case of using optimal K as compared to the case when all users are active, i.e., K = K. 

In the case of high SNR, using G p and (|34l) . the asymptotic expression for the system BER 
can be obtained as 

Pr£r GA W = 0.2g- K T(K + l)(G p Q N )~ K 

(38) 

= 6 ((GpSNR)-*) . 

It follows from (1381) that the asymptotic system BER benefits from the MUDiv power gain 
G~ K = ((K — K)a~ l + l) K at the cost of having the diversity order K < K. 
Finally, inserting (1381) into (|23l) . we obtain that 

d 



lim Kt = argmin 

SNR-+00 K 



arg min \ G P ■ SNR| 



(39) 



dK 
= K + a 

where < a < 1. Since optimal MUDiv K is restricted to be integer, it can be concluded from 
(l39l) that the MUDiv K = K is optimal when SNR— > oo. The latter means that the maximum 
available MUDiv should be used for energy unlimited systems that agrees with known results. 

VI. Extension to multiple antenna MU systems 



The optimization problems proposed in Section [IV] can be extended to multiple antenna MU 
systems, which will also allow to use the benefits of multiple antenna techniques [fT9l , [|20l . 
Toward this end, a generalized expression for the average BER has to be derived. For brevity, 
we consider only the approximate average system BER case. 

Let D denote the multiple antenna diversity order. Then, in the multiple antenna case, the 
degrees of freedom (DOF) of y in (fl4l) extends to 2D, that is, y ~ xId where xId stan ds for 
the Chi-squared distribution with 2D DOF (refer also to Il20l0 . Therefore, for given K and D, 
the expression for pf~ GA (y) in (fl"5l) can be generalized as lfT9l 



-y/Q , u-,J->-i 



P?- bA (y) = a:— 7(Ay/n) A " 1 ^77^ (40) 

where 7(0,6) = t a_1 e'*dt denotes the lower incomplete gamma function [fT8l . 
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Inserting (|40b into (fl"4l) . we obtain the average BER in the multiple antenna case as 

POD 

Pr*- GA (K) = 0.2K / ^-V^+^S'tA^VrOD^dt. (41) 
Jo 

Finally, the optimization problems proposed Section [IV] can be straightforwardly extended to 
the case of multiple antenna MU systems by using (l4Tj) instead of the corresponding BER 
expressions for the single antenna case. As an example, an extension of the problem (1221) to the 
case of multiple antenna MU systems will be investigated numerically in the following section. 

VII. Numerical results 
Consider an MU system with Gray-coded square M-QAM of size M G {4, 64}. Let (Cm,i, Cm,i) 
in m be (1, 1) for M = 4, while {{C M ,i, c M ,i), • • • , (C m ,b, cm, 5 )} be {(7/12, 1/21), (1/2, 3/7), 
(-1/12, 25/21), (1/12, 9 2 /21), (-1/12, 13 2 /21)} for M = 64 OH. Independent log-normal 
distributed shadowing with mean \i — 1 and standard deviation v = 5 is assumed with pathloss 
dB. Considering the approximate BER, i.e., Prf ~ GA (-) = Prf ~ GA (K), the optimal K, i.e., 
K* of (ED), K* dt of d26]>, or K* ds of (EU), are found. The set K = {1, 2, • • • , K}, c = 1, 
a G {1,2,7.81,31.25}, and E d (K)/E f G [a, A" - 1 + a] are used. Note that the parameter 
a = 2 corresponds to the standard case when 280 pilot sub-carriers and 560 data sub-carriers 
are used per one sub-channel in 10 MHz uplink WiMAX (IEEE 802. 16e) (|2T|. For comparisons, 
we also consider other parameter values. For example, the parameter a = 7.81 can be obtained 
by using 32 pilot and 250 data sub-carriers while a = 31.25 results from using 32 pilot and 
1000 data sub-carriers. 

In the case when K = 100 and a = 1, the ratio X(K)/X GA (or, equivalently, E d (K)/Ef A ) 
is set at the values between dB and 20 dB depending on K. The generic GA scheme is also 
depicted for comparison. 

A. Minimizing the system BER 

Example 1: In our first example, we consider the problem (|24|) and the case when for a given 
Vl N G {5, 10} dB and a = 2, total energy grows with K. Note that E^~ GA (or, equivalently, 
the average power Pt(K)) is an increasing function of K. 

Fig. [3] shows Kl of (|24l) versus Pt for various values of Q N . It can be seen from the figure that 
Kl increases with respect to Pt- The latter means that the maximum available MUDiv should 
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be used for energy unlimited systems, while the optimal MUDiv can be significantly smaller 
than the maximum available MUDiv K for energy limited systems. It can also be observed that 
for a given P T and low Q N , the optimal MUDiv K£ is also small and more energy should be 
allocated for actual data transmission E d (K) in order to achieve better BER. Finally, it can be 
also seen in this figure that Kl of (|24l) that minimizes the approximate BER coincides with K£ 
of (1241) that minimizes the exact BER, which validates the use of approximate BER. 

Fig. H] illustrates the impact of Kl on Prjf~ GA (-) versus P T . In this figure, fi^ = 5 dB and 
a = 2 are taken. It can be seen from the figure that Pr^~ GA (-) based on K% is a decreasing 
function of Pt. Moreover, for M = 4 and Pt < 30.5 dB, the generic GA scheme is optimal 
since it provides minimum Prf _GA (-) and K* = K, in this case. However, when P T > 30.5 dB, 
the optimal K% < K is obtained. For example, when M — 4, K% provides 3 dB power gain at 
Pr^ = 10~ 5 as compared to the generic GA. 

Example 2: In the second example, we consider the problem (|24|) and the case when for a 
given maximum achievable MUDiv K, E^~ GA grows with fijv- hi this case, a G {7.8125, 31.25} 
and K = 50 are used. 

Fig. \5\ shows Kl versus Pt. It can be seen from the figure that K\ is an increasing function 
of Pt and it converges to K if more power (energy) is available for all users in the system. The 
convergence rate depends on a and it is higher for larger a and slower for smaller a. Note that 
the practical values of a are smaller than both values tested in this example (see Example 1). 
It can also be observed that for low P T , less K£Ef is required to achieve a better system BER 
than the one achieved if all K users are active. For example, for Pt = 28 dB and a = 7.8125, 
the achieved K£Ef for K£ = 12 is significantly smaller than the one for the generic GA. 

In Fig. |6l the impact of Kl on Pr^~ GA (-) is illustrated versus Pt- A significant power gain 
is provided by the proposed method as compared to the generic GA scheme. For example, in 
the case when a = 7.8125, the use of K£ provides 6 dB power gain at Pr^~ GA (-) = 10~ 4 . A 
significant power gain can be observed even for large a, i.e., a = 31.25. However, regardless 
of a, the aforementioned power gain vanishes and Kl converges to K if Pt — > oo (see also 

Fig. a. 



B. Minimizing the total energy of all users in the system 

Example 3: In the last example, we consider the problems (l26l) and (|27l) for the DT and DS 
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MU systems, correspondingly. The proposed K-GA scheduling scheme based on K* dt of ( |26l ) 
and K* ds of (1271) is compared to the generic GA scheduling scheme. 

Fig. [7J shows the error probability of the proposed K-QK scheduling scheme Prjf~ GA (-) 
averaged over variations of the channel mean u versus A. The parameters K = 100 and BER t = 
IQ~ 3 are taken. The average error probability of the generic GA is computed for two cases with 
and without variations of the channel mean to. In can be seen from the figure that the average 
Pr^~ GA (-) is maintained below the system requirements (i.e., BER t = 10~ 3 ) for the DT MU 
system. For the DS MU system, the average Pr^~ GA (-) is close to the average Pr^f (■) at low 
SNRs since in order to guarantee a given target BER BER 4 , the outage is not allowed even if 
K is not sufficiently large. It can be also seen that as A increases, the DS MU system performs 
closer to the DT MU system. It is because K < K is sufficiently large to guarantee the target 
BER 4 in both cases. 

Based on K* dt and K ds , it can also be seen in Fig. [8] that the average Pt normalized by the 
power required for the generic GA is a decreasing function of A for both the DT and DS MU 
systems. Moreover, the DT MU system requires less power (energy) than the DS MU system 
while satisfying the system requirement on target BER BER t . For example, at A « 10 dB, the 
DS MU systems with BER 4 = 10~ 3 achieves a power saving gain of 10 dB over the generic 
GA for the same average Pr b = 4 x 10~ 3 (see Figs. [7J and [8]). Figs. [7J and [8] also depict that as 
A increases, Pt converges to the power required for the RA scheduling scheme. 

For the multiple antenna case, Fig. [9] shows Pr^~ GA (-) versus the diversity order D for the 
following parameters K = 5, fi^ = 5 dB, a = 1, and M = 4. It is also assumed that the energy 
is distributed according to (|22|) with Prjf~ GA (-) as derived in Section |VT1 It can be seen from 
this figure that our optimal energy distribution gives a boost in the system BER as compared to 
the generic GA scheduling scheme. 

VIII. Conclusions 

A new realistic energy model which describes the distribution of the total finite users' energy 
between scheduling and data transmission stages is developed for the energy limited uplink MU 
wireless systems. MU scheduling algorithms which maximize the MUDiv gain are derived for 
the aforementioned systems to (i) minimize the overall system BER for a fixed total energy of all 
users in the system or (ii) minimize the total energy of all users for fixed BER requirements. It 
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is shown that for a fixed number of available users, an achievable MUDiv gain can be improved 
by activating only a subset of users from the entire set of users. Using asymptotic analysis, it is 
shown that our approach benefits from MUDiv gains higher than that achieved by the generic GA 
algorithm, which is the optimal scheduling method for energy unlimited systems. In particular, 
when minimizing the system BER, it is found that the achieved MUDiv gain is determined by 
(AK — K)a^ + l) -SNR when K is large. Moreover, in the case of high SNR, the MUDiv power 
gain ((K — K)a^ + l) K can be archived while obtaining the diversity order K. Simulation 
results validate our theoretical observations and show that the proposed J^-GA algorithm based 
on optimizing the number of active users provides significant energy gains for energy limited 
MU wireless systems over the generic GA algorithm. 

Appendix: Derivations of (|23l) and ([251) 

Using (fl9l ), the first derivative of Pt^~ (K) in the optimization problem (l24l) can be 
expressed as 

<} Pr£T GA (K) = bB(K, 1 + /(Jf)) + bK^-B (K, 1 + f(K)) . (42) 



dK b ' a y ' v ' J v " dK 
In turn, the first derivative of B (K, 1 + f(K)) with respect to K in (1421) can be written as 

d B(K,1 + f(K)) = A £ t K -\l - t)'Mdt = £ A t K ~\l - t)™dt. (43) 



dK 

or equivalently as 

d 



B(K, 1 + f(K)) = f 1 (1 -t) fiK) t K - l lntdt + f'{K) [\l-tf- l t f(K) \ntdt (44) 
Jo Jo 

where /'(■) denotes the first derivative of /(■) with respect to K and ln(-) stands for the natural 
logarithm.. 

Using the relationship QH (4.253.1)] 

x u ~ l (l - a; r )^ 1 lna; dx = B (u/r, v) {ip(u/r) - ip(u/r + v)}/r 2 (45) 



where ip(z) = ^lnr(2;) denotes the digamma function for z > 0, the first derivative of the beta 
function in (|44|) can be written as 

J^B(K, 1 + f(K)) = B(K, 1 + f(K)) 

x {^(K) - (1 + f{K))^{K + 1 + f{K)) + f{K)^{l + f{K))}. 

8 Note that a logarithm with any basis can replace the natural logarithm. 
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Inserting (|46l) into (1421) . we find that the solution of (|24l) should satisfy the following equation 



B(if,i+/(if)) i+ir{^(K)-(i+f(^))^(^+i+/W)+/'(W(i+/W)} 



0. (47) 



Since in our system model if > 1 and /(if) > 0, it follows from (gZj that B(if, 1 + / (if)) > 0. 
Therefore, the equality B(K, 1 + f(K)) = holds if and only if K goes to infinity. However, 
for K — >■ oo the assumption of the limited total system user energy is violated, and therefore, 
B(K, 1 + f(K)) in (|47l) must always be positive. Thus, the problem of finding the solution of 
(124)) boils down to the problem of finding the number of users which satisfies the following 
equation 

ij(K) — ip(K + 1 + f(K)) + f'(K){i/>(l + f( K )) -HK + 1 + f(K))} + 1/if = 0. (48) 
Using the following expression [fT8l (8.365.3)] 

n-1 

i/j(x + n) =ip(x) + ^(x + iy 1 (49) 
the differences between the digamma functions in (|48l can be represented alternatively as 

K-l 

xP(K+l + f(K)) - ^(1 + f(K)) = £(1 + /(*Q + -1 • (50) 

z=o 
/W 

V(if + 1 + / W) - V(if ) = ^ (if + O" 1 • (51) 

1=0 

Finally, inserting (l50l) and (I5TT) into (|48~1) . the left hand side of (|48l can be rewritten as 

/(*) -, K-l 



z=o z=o ^ v 7 

Therefore, for given and E^~ GA , the optimization problem (|22~1) can be rewritten as 



if* = arg min m(iT) subject to Etf~ = E. (53) 
This completes the derivation. 
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Fig. 1. Total system BERs Pr b , e {K), Pr b , u (K) and Pr b>a (K) for K-GK when K = K = 1, 10, 50. 
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Fig. 5. Optimal number of users Kl versus Pt when K = 50, M = 4, a £ {7.8125, 31.25}. 




Fig. 6. Total system BER using K* b versus P T when if = 50, M = 4, a G {7.8125, 31.25}. 
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Fig. 7. Average Pr 6 of K-GA in DS MU and DT MU systems when fj, = 1, v = 5, BER 4 = 1CT 3 , a = 1, if = 100. 
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Fig. 8. Average P T of K-Gh in DS MU and DT MU systems when (j, = 1, v = 5, K = 100, a = 1, BER t = 10~ 3 . 
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Diversity order, D 



Fig. 9. Impacts of D, i.e., 2D is the DOF resulting from multiple antennas, on Pr^ GA using K% when K = 5, SI at = 5 dB, 
and M = 4. 
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